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I review and discuss a selected sample of recent results in pNRQCD. 

1 Introduction 

Non-relativistic bound-state systems are characterized by, at least, three widely 
separated scales: the mass m of the particle, the (soft) scale associated to its 
relative momentum ~ mv, s<l, and the (ultrasoft) scale associated to its 
kinetic energy ~ mv 2 . In QED and in the perturbative regime of QCD the 
velocity v of the particle in the bound state may be identified with the coupling 
constant. Moreover, the inverse of the size of the system is also of order mv 
and the binding energy of order mv 2 . Indeed, a systematic treatment of non- 
relativistic bound-state systems in the framework of effective field theories 
(EFT), which takes full advantage of the above energy scale hierarchy, was 
initiated in QED 111 and in more recent years remarkable progress has been 
achieved in the analysis of tt threshold productionB. 

For systems made of b and c quarks (I will denote them generically as heavy 
quarkonia: ip, T, B c , ...) non-perturbative contributions may be relevant. By 
comparing the energy level spacings of these systems (see Fig. with the 
heavy-quark masses (e.g. m& ~ 5 GeV and m c ~ 1.6 GeV) we can still argue 
that the data are consistent with a kinetic energy of the bound quark much 
smaller than the heavy-quark mass and, therefore, with a non-relativistic (NR) 
description of the heavy-quark-antiquark system. However, in dependence of 
the specific system, the scale of non-perturbative physics, Aqcd, may turn out 
to be close to some of its dynamical scales. The physical picture, which then 
arises, may be quite different from the perturbative situation. What remains 
guaranteed, also for heavy quarkonia, is that m ^> Aqcd and that at least the 
mass scale can be treated perturbatively, i.e. integrated out from QCD .order 
by order in the coupling constant. The resulting EFT is called NRQCD0. 

A lot of effort has been put over the last two decades in order to find 
the relevant operators, which parameterize the non-perturbative heavy-quark- 
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Figure 1 : The spectra of bb and cc quarkonia normalized with respect to the spin average of 
the Xo(l-P) an d Xc(l-P) states respectively. 

antiquark interaction, once the mass scale has been integrated out. In some 
classical works DOfj'Q these operators were identified with Wilson loop op- 
erators. At the same time, however, the. -celevance of less extended non- 
perturbative objects was pointed out inBtalllJ for situations where the scale 
of non-perturbative physics is of the order mv 2 or smaller. A first non- 
perturbative derivation p£ some heavy quarkonia potentials in the framework 
of NRQCD was done in til. While a full systematic study of the heavy-quark- 
antiquark systems in an EFT framework, which incorporates all the possible 
dynamical situations (at least in pure gluodynamics) and factorizea-the relevant 
non-perturbative operators, has been recently completed inE3H3'0. 

In the following I will discuss the EFT that may be constructed from 
NRQCD by integrating out the scale of the momentum transfer, assumed to 
be the next relevant scale oLihe system. I shall call the obtained EFT, po- 
tential NRQCD (pNRQCD) E3 . In Sec. | I will consider the situation where 
this scale is much bigger than Aqcd, (more specifically I will consider mv 2 
not smaller than Aqcd)- To this situation belong QED bound states (in the 
appropriate gauge-group limit) and what would be tt bound states. In par- 
ticular, I will review the a 5 In a calculation of the quarkonium spectrum and 
some of its implication for the e + e~ — > tt cross section. It is not a priori clear 
to which heavy-quarkonium states these results apply. As a guideline we may 
take the results ofE£l plotted in Fig. |[ Eventually the internal consistency of 
the EFT and the comparison with the experimental data will provide a way 
to discriminate among the different situations. The quarkonium ground-state 
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Figure 2: The size of some heavy quarkonia, as calculated in is shown with respect to 
the quark- antiquark static potential. 

radii, in particular for the T(1S), appear to fall in a region where the potential 
is characterized by a Coulomb-type behaviour. This suggests that a pertur- 
bative treatment at the momentum-transfer scale may be correct. Instead, 
heavy-quarkonium resonances higher than the ground state fall in a region 
where the potential is no longer of the Coulomb-type. This seems to indicate 
that a perturbative treatment of the momentum transfer scale is not allowed 
for them. In Sec. I will consider this last situation. 



2 Quarkonium at the NNNLO 

In this section I shall discuss heavy quarkonium in the dynamical situation 
where rav 2 is not smaller than Aqce>- This means that at the matching scale to 
pNRQCD, vnv > [i > rav 2 , 1 can still assume that (ultrasoft) gluons and quark- 
antiquark states in color-singlet and color-octet configuration exist. What 
would be toponium in tt threshold production and (likely) heavy-quarkonium 
ground states fall in this situation. The aim is to set up the framework for 
an eventual full NNNLO calculation of the heavy quarkonium masses as well 
as the e + e~ — > tt cross section. More explicitly I will give the leading log 
contributions to the NNNLO. 

The pNRQCD Lagrangian in the situation Aqcd < mv 2 , considering only 
the terms relevant to the analysis of the leading-log corrections at the NNNLO 
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Figure 3: The Feynman graphs of the static Wilson loop (left) contributing to the three- 
loop leading logs of the static matching potential and the corresponding graph in pNRQCD 
(right). The double line indicates the octet propagator. 



of the singlet, reads as follows: 
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where r is the relative coordinate, p = — iV r , and S (— S ll c /\/AQ and are 
the singlet and octet field, respectively. All the gauge fields in Eq. (|l|) are 
functions of the centre-of-mass coordinate and the time t only. 

The functions V are the matching coefficients of pNRQCD. Here I am 
interested only in the singlet sector at In a s accuracy. For this purpose 
can be obtained fmm- matching NRQCD to pNRQCD at 0(l/m°) exactly at 
the two-loop levell!3'ta and with leading-log accuracy at the three-loop levelEJ 
(see Fig. ||). The result reads 
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where (3 n are the coefficients of the beta function and the values of a\ and 
0,2 can be found inEj. For the calculation of the matching potential and 
we need to perform the matching.^ exactly at one- loop level En and with 
leading-log accuracy at two-loop levelcj and exactly at tree leveled and with 
leading-log accuracy at one loop leveled respectively. (The leading logs may be 
extracted by considering the ultraviolet divergences of the pNRQCD diagram 



in Fig. ||, when evaluated on the full octet propagator.) The result reads 
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The quarkonium spectrum at leading log accuracy of the NNNLO, in the 
situation Aqcd < mi) 2 , is given by EH 
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where E n — — to f „ s , /i G = p 2 /to + V'J '' and the states In, i) are the eigen- 

states of the Hamiltonian p 2 /m + V^°\ The [i dependence of the first line of 
Eq. cancels against the ultrasoft contributions of the second line, which 
corresponds to the pNRQCD diagram of Fig. || when evaluated on the full 
octet propagator i/(E — h Q ). 

In the situation where Aqcd < tow 2 the correlator (E a (t)4>(t, 0)^E b (0)) 
can be calculated perturbatively (a one- loop calculation of it is inEJ). An 



explicit expression of Eq. (jll]) at order In a s is given in EI In Ej also the 
ultrasoft corrections to the wave- functions in the origin have been calculated: 

C F C A + \c 2 A \\n 2 — . 

3 J a s 

(12) 

From these the leading-log correction to the NNNLO of the e + e — ► tt cross 
section has been calculated. In Fig. || R(E) = a(e + e~ — ► ti)/a(e + e~ — > /i + /i~) 
is shown at NNLO and with the leading log correction included. The aim of 
such an analysis, once the complete NNNLO will be calculated, is to reach 
a 50 MeV sensitivity on the top quark mass from^he t-t cross-section near 
threshold to be measured at a Next Linear Collider □. 
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Figure 4: R(E) versus E for the parametpn choice rat = 175 GeV, Tt = 1A3 GeV and 
a s (Mz) = 0.118. The plot is taken from E3. The dotted line corresponds to the NNLO 
calculation, the solid line includes the leading log NNNLO. 



For the T(IS') the order <x! lna s correction evaluated from Eq. (^TJ) reads 

1730 

8E 10A = ——m b a*(r)a s ln(l/a s ). (13) 

S17T 

This correction has been considered inEl'0. In dependence of the (ultrasoft) 
scale at which a s is calculated it may be as large as 80—100 MeV. It is not clear, 
up to now, if the size of this correction should be taken as a serious estimate 
of the complete order a^, or if it will largely cancel against the remaining q s 5 
contributions (still unknown). It is also possible that this is a signal of the 
renormalon of order Aq CD r 2 affecting thematic potential, which has proved 
to be related to this kind of corrections int3. 



Despite some remarkable progress achieved recently in increasing our know- 
ledge aLperturbative corrections either by resumming potentially large loga- 
rithms E3 or by consideringJn the bottomonium system the effects due to the 
finite charm quark mass E3, the real challenge for heavy quarkonia remain 
non-perturbative contributions. The uncertainty related to them is usually 
believed to be of 100 MeV for the T(15) and of several hundreds MeV for 
the J/if) and, therefore, it dominates over higher perturbative corrections, 
once a renormalon free mass definition has been used. The leading non- 
perturbative contributions to the spectrum can be also read from Eq. (pT|). 
For the general case Aqcd < mv 2 they are encoded into the gluonic correlator 

(E a (t)4>(t, 0)^E b (0)), which may be expanded in terms of lncal condensates 
in the situation Aqcd <C mv 2 . A discussion can be found inE3. It is not clear 
what situation applies to the physical systems of interest. If the ground states 
of bottomonium and charmonium fall into the situation Aqcd < mv 2 , which 
is likely, a study of these systems using Eq. ( |TT| ) and one of the_uarameteriza- 
tion of the gluonic correlator suggested by sum-rule palculationsEj, by different 
lattice simulations EjEj or by QCD vacuum modelst^l is timely. 
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Figure 5: Quarkonium-quarkonium scattering and extraction of the van der Waals potential. 



2.1 Quarkonium-quarkonium scattering 

The above EFT approach may be further pursued if the heavy quarkonium sys- 
tem interacts with other systems so that scales smaller than the binding energy, 
E blnd , are present. This is the situation that may happen in the scattering of 
heavy-quarkonium states if the energy of the hadron, £? had , is much smaller 
than E bmd o. By integrating out from the scattering amplitude the higher 
energy scales, we may get a suitable definition of the quarkonium-quarkonium 
potential (see Fig. ||). The situation is similar to the matching to pNRQCD 
discussed above. More specifically the quarkonium-quarkonium van der Waals 



potential is given by 
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being R the relative coordinate of the quarkonia. The coefficients a 2 - 7 come 
from the matching to pNRQCD: 
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Notice that in this case the relevant non-perturbative operator is a four- 
chromoelectric-field correlator. Applications to the J/tp — J/tp scattering have 
been recently discussed inEj. 

3 Long-range quarkonium 

For higher heavy-quarkonium states the Coulombic Bohr radius tends to be- 
come large and the perturbative matching, which led to pNRQCD in the above 
formulation, is no longer justified^] However, the success of traditional poten- 
tial models seems to suggest, that a NR description of these systems may still 
hold (for some reviews seeEj'EZl). Therefore, one may still think to follow the 
same procedure discussed in the previous section for the perturbative case and 
integrate out the scale of the momentum transfer in order to get what would be 
pNRQCD in this situation. The result is a NR quantum-mechanical descrip- 
tion of heavy quarkonium fully derived from (NR)QCD via a non-perturbative 
matching. I will assume that the matching between NRQCD and pNRQCD 
can be performed order by order in a 1/m expansion. While this can be jus- 
tified within a perturbative framework, in a non-perturbative situation, one 
may question on its validity. For instance, a case where certain degrees of 
freedom cannot be integrated out in this way has been considered inE£l. This 
point surely deserves further studies. It is also relevant in order to establish 
the proper power counting of NRQCD, on which I will comment later on. 

In order to identify the degrees of freedom of pNRQCD in the situation 
where the momentum transfer of the system is close to Aqcd , I first consider 

the NRQCD Hamiltonian up to order 1/m 2 , H = H (0) + where 



n=l,2 



H (0) = J rf 3 xi (LFIT + B a B a ), II a is the canonical momentum conjugated 

b It is debated if for the charmonium ground state the scale of the momentum transfer may 
or not be integrated out perturbatively. Here, I note that, if this cannot be done, also the 
mentioned approach to the J/rjj — J/tp scattering needs to be revised. 



to the gauge field A £ j_and the terms iJW and may be read off from the 
NRQCD LagrangianEil In the static limit the one-quark-one-antiquark sector 
of the Fock space may be spanned by |n; xi, x 2 )(°) = ip^ (xi)x£( x 2)|™; x i, x 2)^°\ 
where |n;xi,x 2 )(°) is a gauge-invariant eigenstate (up to a phase) of 
with energy E n (xi,X2), and Xc( x ) = *c 2 x*(x). \n; xi, x 2 )(°) encodes the 
gluonic content of the state, i.e. it is annihilated by Xc( x ) and - 0( x ) (Vx). The 
positions Xi and x 2 of the quark and antiquark respectively are good quantum 
numbers for the static solution |n;xi,x 2 )(°) (but will be used also to label 
the eigenstates of H); n generically denotes the remaining quantum numbers, 
which are classified by the irreducible representations of the symmetry group 
Dooh (substituting the parity generator by CP). I also choose |n;xi,x 2 )(°) 
to be invariant under time inversion. The ground-state energy Eq (xi,X2) 
can be associated (in some specific situation) to the static potential of the 
heavy quarkonium. The remaining energies En (xi,x 2 ), n ^ 0, are usually 
associated to the potentials describing heavy hybrids or heavy quarkonium (or 
other heavy hybrids) plus glueballs. They can be computed on the lattice (see, 
for instance, Eac3). Beyond the static limit, but still working order by order in 
1/m, the eigenvalues £' n (xi,x 2 ;p) of the Hamiltonian H, up to 0(l/m 2 ), are 
given by 
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Eq corresponds to the quantum-mechanical Hamiltonian of the heavy quarko- 
nium (in some specific situation). The other energies E n , for n>0, are related 
to the quantum-mechanical Hamiltonian of higher gluonic excitations between 
heavy quarks. Explicit expressions for the energies E n , obtained from the 
above formula, can be found inE3'cj. 

Let me now assume that, because of a mass gap in QCD, the energy split- 
ting between the ground state and the first gluonic excitation is larger than 
mv 2 (see also the data reported by G. Bali at this conference), and, because 



of chiral symmetry breaking of QCD, Goldstone bosons (pions/kaons) appear. 
Hence, in this situation, the states with ultrasoft energies (i.e. the degrees 
of freedom of pNRQCD) would be the ultrasoft excitations about the static 
ground state, which we call the singlet, plus the Goldstone bosons. If one 
switches off the light fermions (pure gluodynamics) , only the singlet survives 
and pNRQCD reduces to a pure two-particle NR quantum-mechanical system. 
Therefore, under the assumption of the validity of the 1/m expansion (in the 
matching) and of the existence of a mass gap between the singlet and the other 
gluonic excitations between heavy quarks, we obtain the typical situation de- 
scribed by potential models. More specifically, in terms of the NRQCD states 
discussed above, this means that only |0; Xi, x 2 )(°) is kept as an explicit degree 
of freedom, whereas |n; Xi, x 2 )(°) with n ^ are integrated out. |0;xi,x 2 )(°) 
provides the only dynamical degree of freedom of the theory. It is described 
by means of a bilinear colour singlet field, S^xi, x 2 , i), with the same quan- 
tum numbers and transformation properties under symmetries. In the above 
situation, the Lagrangian of pNRQCD reads 



-CpNRQCD = >S T I «? _ h s {xi,x 2 ,p)jS, (16) 

p 2 p 4 y(l) y{2) 

where h s = ? + ' H 1 5- is the Hamiltonian of the singlet 
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and may be identified through the matching condition 

£b(xi,x 2j p) = /i s (xi,x 2 ,p). (17) 

I note that, if other ultrasoft degrees of freedom, apart from the singlet, exist, 
they may be added systematically to the above Lagrangian in an analogous 
way as done in the perturbative situation (see Eq. ([!])). For what concerns 
the effects on the computation of the potentials, since we are integrating over 
all the states, in the situation where some of them, different from the singlet, 
are ultrasoft, wc would just need to subtract their contribution later on. 

InEj'Eil the matching of NRQCD to pNRQCD has been performed up to 
order 1 /to 2 and the above potentials have been obtained explicitly in terms of 
Wilson loops □. For the static potential the result reads 

V (0) = lim Iwwb), (18) 

where Wn is a rectangular WilsoxLlopp of dimension rxT. The 1 /to and 1 / m 2 
potentials may be read off fromlija after the identifications (cf. Eq. (|J)): 



Vg=Vg?\ Vg» = Vg\ (19) 

Having expressed the non-perturbative dynamics of the heavy-quark potentials 
in terms of Wilson loops is extremely convenient, for these quantities may be 
calculated directly in lattice simulations. tiL Moreover, these operators can be 
also calculated in QCD vacuum modelse3'E3, providing a way to check, directly 
on the phenomenology, assumptions on the structure of the QCD vacuum. In 
particular, it would be of interest to see what the vortices picture of the QCD 
vacuum, which has received so much attention in recent yearsE-3, may predict on 
these correlators. Finally, I would like to stress that the obtained expressions 
for the potentials are also correct perturbatively at any order in a s . 



3.1 The NRQCD power counting 

An important issue, once the EFT Lagrangian has been calculated through 
the matching procedure, is to establish its power counting in order to calculate 
physical observables.| Establishing the power counting of pNRQCD in the 
non-perturbative regime is, however, not only important by itself. It may also 
serve to establish the non-perturbative power counting of NRQCD, which is an 
important source of informatiorLon the spectrum of excited quarkonium states 
(but also heavy-light mesons) l3 The power counting usually adopted there 
and discussed, for instance, inEJ is inherited from the perturbative regime. 
However, there is no certainty that this is the suitable one for calculating 
higher quarkonium states, since in the non-perturbative regime different power 
countings are, in principle, possible. The above formulation of pNRQCD has 
translated the problem of the NRQCD power counting to obtaining the power 
counting of the different potentials. This is expected to be of some advantage: 
1) because the power counting of pNRQCD is simpler and quantum-mechanical 
arguments, like the virial theorem, may be more properly formulated in this 
context; 2) because all the potentials are expressed in terms of Wilson loops, 
for which there arc or there will be direct lattice measurements. 

As an example of power counting in pNRQCD we can assume that the 
potentials scale with mv. By definition the kinetic energy counts as mv 2 . 
would count as mv, if the virial theorem would not constrain it to count also 
as mv 2 . In the perturbative case this extra v suppression comes from the 
factor a s ~ v in the potential (see Eq. (§)). From our general assumption 
V^/m scales like mv 2 . Therefore, it could be in principle as large as V^°\ 



c It should be stressed that the way the matching procedure is organized, e.g. as an l/m 
expansion, may be different from the power counting of the EFT. 



This makes a (still to do) lattice evaluation of this potential quite interesting. 
Perturbatively, due to the factor (see Eq. (g)), it is 0(mv 4 ). For what 
concerns the 1/m 2 potentials, they are in this scheme of order mv 3 . However, 
also here several constraints apply. Terms involving V(°) are suppressedJiv 
an extra factor v, due to the virial theorem. General symmetry relations ErIZI, 
also further constrain the power counting. Finally, some of the potentials are 
0(a s ) suppressed in the matching coefficients inherited from NRQCD. 

Acknowledgments 

I thank Nora Brambilla for reading the manuscript and comments. I thank the 
organizers for the invitation, the Alexander von Humboldt foundation for sup- 
port and the CFIF at the Instituto Superior Tecnico of Lisbon for hospitality 
during the last stage of this work. 

References 

1. W. E. Caswell and G. P. Lepage, Phys. Lett. B167, 437 (1986). 

2. A. Hoang at al., Eur. Phys. J. direct C3, 1 (2000). 

3. B. A. Thacker and G. P. Lepage, Phys. Rev. D43, 196 (1991). 

4. K. G. Wilson, Phys. Rev. D10, 2445 (1974). 

5. L. Susskind, Les Houches lectures (1976); W. Fischler, Nucl. Phys. 
B129, 157 (1977); L. S. Brown and W. I. Weisberger, Phys. Rev. D20, 
3239 (1979). 

6. E. Eichten and F. L. Feinberg, Phys. Rev. D23, 2724 (1981); M. E. 
Peskin, in Proceeding of the 11th SLAC Institute, SLAC Report No. 
207, 151, edited by P. Mc Donough (1983); D. Gromes, Z. Phys. C26, 
401 (1984); in Spectroscopy of light and heavy quarks, ed. by U. Gastaldi, 
R. Klapisch and F. Close, (Plenum Press, New York, 1989). 

7. A. Barchielli et al., Nucl. Phys. B296, 625 (1988); (E) 303, 752 (1988); 
A. Barchielli et al., Nuovo Cimento 103, 59 (1990); N. Brambilla et al., 
Phys. Rev. D50, 5878 (1994). 

8. M. B. Voloshin, Nucl. Phys. B154, 365 (1979); Sov. J. Nucl. Phys. 36, 
143 (1982); H. Leutwyler, Phys. Lett. B98, 447 (1981). 

9. D. Gromes, Phys. Lett. B115, 482 (1982). 

10. I. I. Balitsky, Nucl. Phys. B254, 166 (1985). 

11. Y. Chen et al., Phys. Rev. D52, 264 (1995); N. Brambilla and A. Vairo, 
Nucl. Phys. B(Proc. Suppl.) 74, 201 (1999). 

12. N. Brambilla et al., Phys. Rev. D60, 091502 (1999); Nucl. Phys. B566, 
275 (2000). 

13. N. Brambilla et al., |hcp-ph/0002250| . 



14. A. Pineda and A. Vairo, |hcp-ph/0009145| . 

15. S. Godfrey and N. Isgur, Phys. Rev. D32, 189 (1985). 

16. A. Pineda and J. Soto, Nucl. Phys. B(Proc. Suppl.) 64, 428 (1998). 

17. A. Billoire, Phys. Lett. B92, 343 (1980). 

18. Y. Schroder, Phys. Lett. B447, 321 (1999); M. Peter, Phys. Rev. 
Lett.78, 602 (1997). 

19. A. Duncan, Phys. Rev. D13, 2866 (1976). 

20. N. Brambilla et al., Phys. Lett. B470, 215 (1999). 

21. S. N. Gupta et al., Phys. Rev. D26, 3305 (1982); S. Titard and F. J. 
Yndurain, Phys. Rev. D49, 6007 (1994). 

22. B. Knichl and A. Penin, Nucl. Phys. B563, 200 (1999). 

23. B. Kniehl and A. Penin, Nucl. Phys. B577, 197 (2000). 

24. M. Eidemul ler and M. Jamin , Phys. Lett. B416, 415 (1998). 

25. M. Beneke, |hep-ph/9911490|. 

26. F. Yndurain~ plcp-ph/0007333| ; |icp-ph/0002237| . 

27. M. Luke et al., Phys. Rev. D61, 07402 5 (2000); A. Ma nohar and I. 
Stewart, Phys. Rev. D62, 014033 (2000); |hep-ph/0003107| ; A. Manohar 
et al, Phys. Lett. B486, 400 (2000) 



J. Soto and A. Pineda, bep- 



ph/0007197. 



28. 
29. 
30. 
31. 
32. 
33. 

34. 



35. 
36. 
37. 

38. 
39. 

40. 
41. 
42. 



A. Hoang, |hcp-ph/0008T02| ; D. Eiras and J Soto, |icp-ph/0005066 . 

N. Brambilla and A. Vairo, |hcp-ph/0004192| . 

H. G. Dosch et al., Phys. Lett. B452, 379 (1999). 

M. D'Elia et al., Phys. Lett. B408, 315 (1997). 

G. S. Bali et al., Phys. Lett. B421, 265 (1998). 

M. Baker et al., Phys. Rev. D58, 034010 (1998); D. Antonov, JEEP 07, 
055 (2000). 

M. E. Peskin, Nucl. Phys. B156, 365 (1979); G. Bhanot and M. E. Pe- 
skin, Nucl. Phys. B156, 391 (1979); see also B. R. Holstein contribution 
at INT 00-2 (2000). 

H. Fujii and D. Kharzeev, Phys. Rev. D60 1 14039 (1999). 
N. Brambilla an d A. Vairo, |hcp-ph/990433C . 

F. J. Yndurain, |hcp-ph/9910399| ; Nucl. Phys. B(Proc. Suppl.) 64, 433 
(1998). 

D. Eiras an d J. Soto, Phys. R ev. D61, 114027 (2000). 

C. Michael, hep-ph/9809211 and at this conference; K. J. Juge et al., 



hep-lat/9809015 



N. Brambilla, Nucl. Phys. B(Proc. Suppl.) 86, 389 (2000) 
G. S. Bali et al., Phys. Rev. D56, 2566 (1997). 
N. Brambilla and A. Vairo, Phys. Rev. D55, 3974 (1997); 
al., Phys. Lett. B389, 577 (1996); L. Fulcher, Phys. Rev. 



M. Baker et 
D62, 094505 



(2000) and at this conference. 

43. N. Brambilla, jhep-ph/9809263 . 

44. H. Rcinhardt at this conference; S. Olejnik at this conference. 

45. H. Shanahan at this conference; S. Collins at this conference. 

46. G. P. Lepage et al., Phys. Rev. D46, 4052 (1992). 



